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Three-dimensional  topological  Weyl  semimetals  can  generally  support  a  zero-dimensional  Weyl 
point  characterized  by  a  quantized  Chern  number  or  a  one-dimensional  Weyl  nodal  ring  (or  line) 
characterized  by  a  quantized  Berry  phase  in  the  momentum  space.  Here,  in  a  dissipative  system 
with  particle  gain  and  loss,  we  discover  a  new  type  of  topological  ring,  dubbed  Weyl  exceptional 
ring  consisting  of  exceptional  points  at  which  two  eigenstates  coalesce.  Such  a  Weyl  exceptional 
ring  is  characterized  by  both  a  quantized  Chern  number  and  a  quantized  Berry  phase,  which  are 
defined  via  the  Riemann  surface.  We  propose  an  experimental  scheme  to  realize  and  measure  the 
Weyl  exceptional  ring  in  a  dissipative  cold  atomic  gas  trapped  in  an  optical  lattice. 


Recently,  condensed  matter  systems  have  proven  to  be 
a  powerful  platform  to  study  low  energy  gapless  particles 
by  using  momentum  space  band  structures  to  mimic  the 
energy-momentum  relation  of  relativistic  particles  [1,  2] 
and  beyond  [3-6].  One  celebrated  example  in  three  di¬ 
mensions  is  the  zero-dimensional  Weyl  point  [7-18]  de¬ 
scribed  by  the  Weyl  Hamiltonian,  which  has  been  long 
sought-after  in  particle  physics  but  only  experimentally 
observed  in  condensed  matter  materials  [19-21].  Such  a 
Weyl  point  can  be  viewed  as  a  magnetic  monopole  [22] 
in  the  momentum  space  and  possesses  a  quantized  Chern 
number  on  a  surface  enclosing  the  point.  Another  exam¬ 
ple  is  the  one-dimensional  Weyl  nodal  ring  [3,  23-25], 
which  has  no  counterpart  in  particle  physics.  It  can  be 
regarded  as  the  generalization  of  zero-dimensional  Dirac 
cones  in  two-dimensional  systems,  such  as  in  graphene, 
to  three-dimensional  systems.  Such  a  nodal  ring  has  a 
quantized  Berry  phase  over  a  closed  path  encircling  it 
but  does  not  possess  a  nonzero  quantized  Chern  number. 
This  leads  to  a  natural  question  of  whether  there  exists  a 
topological  ring  exhibiting  both  a  quantized  Chern  num¬ 
ber  and  a  quantized  Berry  phase  in  the  momentum  space. 

So  far,  studies  on  those  gapless  states  focus  on  closed 
and  lossless  systems.  However,  particle  gain  and  loss  are 
generally  present  in  natural  systems.  Such  systems  can 
often  be  described  by  non-Hermitian  Hamiltonians  [26— 
29],  which  are  widely  applied  to  many  different  sys¬ 
tems  [30-40].  Due  to  the  non-Hermiticity,  eigenvalues 
of  the  Hamiltonian  are  generically  complex  unless  the 
Vl~  symmetry  [41]  is  conserved  and  the  imaginary  part 
of  energy  is  associated  with  either  decay  or  growth.  An¬ 
other  intriguing  feature  of  a  non-Hermitian  system  is  the 
existence  of  exceptional  points  (EPs)  [26-29]  at  which 
two  eigenstates  coalesce  and  the  Hamiltonian  becomes 
defective,  leading  to  many  novel  phenomena,  such  as  loss- 
induced  transparency  [30],  single-mode  lasers  [36,  37], 
and  reversed  pump  dependence  of  lasers  [33] . 

In  this  paper,  we  investigate  a  system  of  Weyl  points 
in  the  presence  of  a  spin-dependent  non-Hermitian  term 
and  find  a  Weyl  exceptional  ring  composed  of  EPs.  In 
stark  contrast  to  a  Weyl  nodal  ring  [3,  24,  25],  which 
does  not  have  a  nonzero  Chern  number,  remarkably,  this 


ring  exhibits  a  nonzero  quantized  Chern  number  as  long 
as  the  integral  of  the  Berry  curvature  is  evaluated  over 
a  surface  (labeled  by  S)  that  encloses  the  whole  ring. 
Since  energy  is  multi-valued  in  the  complex  parameter 
space  due  to  its  square  root  form,  a  state  on  the  surface 
S  may  be  defined  over  the  Riemann  surface,  on  which  a 
function  is  single  valued.  On  the  other  hand,  the  Chern 
number  is  zero  when  the  surface  S  does  not  enclose  any 
part  of  the  ring  even  when  it  is  located  inside  it.  Besides 
the  Chern  number,  such  a  Weyl  exceptional  ring  has  a 
quantized  Berry  phase  over  a  trajectory  encircling  the 
ring  twice,  instead  of  once  in  the  case  of  the  Weyl  nodal 
ring.  Furthermore,  we  propose  a  feasible  scheme  to  engi¬ 
neer  and  probe  the  Weyl  exceptional  ring  in  a  dissipative 
ultracold  atomic  gas.  In  such  a  system,  we  find  that  the 
Fermi  arc  can  still  exist  but  is  suppressed,  even  though 
the  Weyl  point  transforms  into  a  ring. 

Toy  model  of  Weyl  exceptional  ring — Near  a  Weyl 
point  in  the  momentum  space,  a  system  can  be  described 
by  the  Weyl  Hamiltonian  Hw  =  ±  vvkuav,  where  av 
represent  Pauli  matrices  and  ±  the  chirality.  For  clar¬ 
ity,  we  consider  the  positive  chirality  and  choose  vv  =  1 
hereafter.  In  the  presence  of  a  non-Hermitian  term  i^az 
(7  >  0)  associated  with  particle  gain  for  spin  up  atoms 
and  loss  for  spin  down  ones,  the  Hamiltonian  becomes 

H(k)=  kvav  +  *707,  (1) 

v=x,y,z 

taking  the  energy  unit  to  be  1.  The  eigenvalues 
are  Eg( k)  =  i/fc2  —  72  +  2zfczy  =  \/ A(k)elB/2,  where 
A(k)  =  yj (fc2  —  72)2  +  4fc2y2  with  k2  =  fc2  +  fc2  +  fc2, 
and  0  is  defined  via  cosd  =  (fc2  —  72)/A(k)  and  sind  = 
2fcz7/A(k).  Here,  6  is  used  to  label  two  branches,  given 
that  eifl/2  gains  a  minus  sign  upon  0  — ►  6  +  2i r,  corre¬ 
sponding  to  the  other  band.  In  the  absence  of  7,  energy 
of  both  bands  is  real  and  two  bands  touch  at  k  =  0 
with  linear  dispersion  along  all  three  momentum  direc¬ 
tions.  In  this  case,  0  takes  only  two  nonequivalent  dis¬ 
crete  values:  0  and  27r  (corresponding  to  two  distinct  and 
separate  bands).  When  7  >  0,  the  eigenvalues  become 
complex,  and  the  single  touching  point  morphs  into  a 
Weyl  exceptional  ring  in  the  kz  =  0  plane  characterized 
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FIG.  1.  (Color  online)  Energy  spectra  and  the  Riemann  surface  of  the  toy  model  in  Eq.  (1).  Spectra  with  respect  to  kx  and  ky 
for  k,  =  0  in  (a)  (real  parts)  and  (b)  (imaginary  parts).  Real  (c)  and  imaginary  parts  (d)  of  the  Riemann  surface  as  a  function 
of  ky  and  kz  for  kx  =  0.  In  (c)  and  (d),  the  color  represents  the  strength  of  9  mod  4-7T,  and  the  red  tube  arrow  shows  a  path 
from  6  =  0  to  9  =  47i\ 


by  fc2  +  ky  =  y2.  On  this  ring,  both  the  real  and  imagi¬ 
nary  parts  of  the  eigenvalues  vanish  [shown  in  Fig.  1(a) 
and  (b)]  and  two  eigenstates  coalesce  into  a  single  one 
(different  from  the  case  of  degeneracy).  Fig.  1(a)  and  (b) 
also  illustrate  that  in  the  kz  =  0  plane  energy  is  purely 
real  outside  the  ring  (with  conserved  VT  symmetry)  and 
purely  imaginary  inside  it  for  this  simple  model.  Inter¬ 
estingly,  9  takes  continuous  values  from  0  to  47t  (9  and 
9  +  4-7T  are  equivalent)  and  gains  27t  when  a  state  trav¬ 
els  through  the  ring  and  returns,  ending  up  at  another 
state  with  opposite  energy,  arising  from  the  role  of  branch 
points  that  the  Weyl  exceptional  ring  plays. 

In  complex  analysis,  besides  using  branch  cuts,  an  al¬ 
ternative  visual  representation  to  depict  a  multi-valued 
function  is  the  Riemann  surface,  a  two-dimensional  (2D) 
manifold  that  wraps  around  the  complex  plane  infinite 
(noncompact)  or  finite  (compact)  number  of  times.  Be¬ 
fore  we  discuss  the  topology  of  the  Weyl  exceptional  ring, 
let  us  first  focus  on  the  definition  of  a  closed  2D  surface 
S  via  the  Riemann  surface.  In  Fig.  1(c)  and  (d),  we  plot 
the  Riemann  surface  of  Eg  for  kx  =  0  (the  color  repre¬ 
sents  the  strength  of  9  mod  47t),  showing  that  energy  is 
single-valued  on  the  surface,  which  connects  the  different 
bands.  Given  the  single  value  property,  we  define  each 
state  on  S  to  be  living  on  the  Riemann  surface.  For  ex¬ 
ample,  if  we  consider  a  state  at  ko  with  9q,  any  other 
states  on  the  surface  S  can  be  obtained  by  starting  from 
this  state  and  travelling  on  the  momentum  space  surface 
S  while  keeping  Eg{ k)  on  the  Riemann  surface. 

With  the  proper  definition  of  a  closed  2D  surface,  we 
can  characterize  the  topology  of  a  Weyl  exceptional  ring 
by  the  Chern  number  on  the  surface  based  on  two  ap¬ 
proaches:  the  integral  of  spin  vector  fields  and  the  Berry 
curvature.  For  the  former,  the  Chern  number  is  given 
by  [42] 

Nl  =  h ld' x  a£) dUldU2-  (2) 

which  characterizes  the  number  of  times  that  the  spin 
field  d g  =  Ei/=iSz(<Ti')e>'  wraps  around  a  closed  sur¬ 


face  S  parametrized  by  (zii,^).  Here,  e„  denotes 
the  unit  vector  along  the  v  direction  and  (ay)  = 
(ug(k)\au\ug(k))  with  \ug(k))  being  the  normalized  right 
eigenstate  of  H( k)  [i.e.,  H(k)\ug(k))  =  Eg(k)\ug{k)) 
and  (ug(k)\ug(k)}  =  1].  Direct  calculations  show  that 
N:i  =  ±1  when  the  surface  S  encloses  the  whole  ring  as 
shown  in  Fig.  2(a),  while  N3  =  0  when  it  does  not  enclose 
any  part  of  the  ring  [shown  in  Fig.  2(b)]. 

Analogous  to  the  scenario  without  decay  [43],  we  may 
also  define  the  first  Chern  number  via  the  Berry  curva¬ 
ture 

C2  =  - j>sn°(k)-dS,  (3) 

where  Clg(k)  =  i(Vku,g(k) |  x  \X7kug(k)}  is  the  Berry  cur¬ 
vature.  Our  calculations  show  that  C3  =  ±1  when  the 
surface  S  encloses  the  Weyl  exceptional  ring  and  C2  =  0 
otherwise,  suggesting  that  the  topological  charge  is  en¬ 
tirely  carried  by  the  ring. 

The  physical  meaning  of  the  Berry  curvature  in  this 
system  can  be  understood  from  the  following  semiclassi- 
cal  equation  under  an  external  gradient  force  F  (see  the 
supplemental  material  for  derivation) 

rc  =  dkcE(kc)  -  kc  x  S 1g(kc),  (4) 

hkc  =  F,  (5) 

where  E(kc)  =  Re[F;0(kc)]  +  A0(kc)  •  kc,  Ae(kc)  = 
Re[A0(kc)  —  A0(kc)]  with  the  Berry  connection  be¬ 
ing  A0(k)  =  *(u0(k)|dku0(k))  and  A0(k)  = 

*(u0(k)|<9kW0(k)),  where  (ug(k) |  is  the  normalized  left 
eigenstate  of  H  [i.e.,  ( ug(k)\H(k )  =  (ug(k)\Eg(k)  and 
(n0(k)|u0(k))  =  1];  rc  and  kc  are  the  center  coordinate 
of  a  wave  packet  in  the  real  space  and  momentum  space, 
respectively.  Clearly,  the  Berry  curvature  plays  the  same 
role  as  in  the  traditional  semiclassical  equation  in  a  closed 
system  [44],  However,  in  this  open  system,  the  equation 
includes  a  term  that  effectively  modifies  the  energy  spec¬ 
tra,  resulting  from  the  difference  between  left  and  right 
eigenstates,  a  feature  in  a  non-Hermitian  Hamiltonian. 
Without  F,  the  group  velocity  is  dictated  by  the  real 
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part  of  the  spectra,  which  implies  that  inside  the  Weyl 
exceptional  ring  in  the  kz  =  0  plane,  the  group  velocity 
vanishes. 

Other  than  the  Chern  number  on  the  surface,  there 
also  exists  a  quantized  Berry  phase  characterizing  the 
Weyl  exceptional  ring,  defined  as 


Ci=<£  z({tfl(k)|5kU0(k))  •  dk,  (6) 

J2C 


where  the  path  2 C  travels  across  the  ring  twice  along  the 
Riemann  surface  so  that  the  state  returns  to  the  origi¬ 
nal  one  after  the  entire  trajectory  as  shown  in  Fig.  1(c) 
and  (d).  Direct  calculations  yield  C\  =  ±7r,  consistent 
with  the  result  for  a  single  EP  [26].  This  Berry  phase 
is  different  from  that  of  a  Weyl  nodal  ring  in  which  the 
quantized  Berry  phase  is  obtained  when  the  trajectory 
encircles  the  ring  once  [3,  24,  25]. 

Realization  in  dissipative  cold  atomic  gases — To  realize 
the  Weyl  exceptional  ring  in  cold  atoms,  we  consider  the 
following  model 


H 


[(^+*7)4,  ,^zCkz,x+  (-^4„x2fc..: 


x+ae^ 


kz  ,x 


u=x,y 


+  (-iy*+lv  JSOvc\  xC>VCfe*,x+ae„  +  H.C.)  +  h0  ,  (7) 


where  x  =  jxaex  +  jyaey  (with  a  being  the  lat¬ 
tice  constant)  labels  the  location  of  sites,  c[,  x  = 

(  4.,x,t  4„x4  )  with  4*,x,<7  (cfea>x,a)  being  the  creation 
(annihilation)  operator,  J  and  Jsov  ( JsOx  =  —JsOy  = 
J so)  stand  for  the  tunneling  and  spin-orbit  coupling 
strength,  ho  =  [~i"i+h2k1/ {2m)]c'k^  xCkz,x,  withy  denot¬ 
ing  the  decay  strength,  and  hz  =  akz  +  hz  is  the  effective 
Zeeman  field  with  a  =  h2ki,z/  (2m)  where  k^z  depends 
on  the  wave  vector  of  Raman  laser  beams  along  the  2 
direction,  m  is  the  mass  of  atoms,  and  hz  the  Zeeman 
field  proportional  to  the  two-photon  detuning.  Here,  we 
consider  the  atoms  to  be  trapped  in  an  optical  lattice  in 
the  x  and  y  directions  while  there  is  no  lattice  along  the 
z  direction. 

Without  7,  this  Hamiltonian,  which  has  two  Weyl 
points  and  a  fourfold  degenerate  point,  can  be  experimen¬ 
tally  engineered  by  coupling  two  hyperfine  states  with 
two  pairs  of  Raman  laser  beams  in  cold  atom  optical  lat¬ 
tices  [45].  To  generate  the  decay  term  representing  an 
atom  loss  —  2iy  for  spin  down  atoms,  one  may  consider 
using  a  resonant  optical  beam  to  kick  the  atoms  in  the  |j,) 
state  out  of  a  weak  trap  as  shown  in  Fig.  2(d),  which  has 
been  experimentally  realized  in  6Li  [40].  Alternatively, 
one  may  consider  applying  a  radio  frequency  pulse  to  ex¬ 
cite  atoms  in  the  |j_)  state  to  another  irrelevant  state  |3), 
leading  to  an  effective  decay  for  spin  down  atoms  when 
atoms  in  |3)  experience  a  loss  by  applying  an  anti-trap. 

To  see  the  energy  spectra,  we  write  down  the  Hamil¬ 
tonian  in  the  momentum  space, 


H{ k)  =  (hz  +  i^)az  -  htrx  +  Ty(-bxax  +  byay),  (8) 


FIG.  2.  (Color  online)  (a)  A  surface  enclosing  a  Weyl  ex¬ 
ceptional  ring  and  (b)  a  surface  located  inside  the  ring,  (c) 
Lattice  structure  in  the  (x,  y)  plane,  (d)  Schematic  of  trapped 
atoms  being  kicked  out  by  a  resonant  optical  beam  (denoted 
by  the  grey  arrow). 


in  the  basis  *F(k)T  with  H7  (k)  = 

(  eifc*aikT  eik*aAkl  HkT  Bkl  ),  where  Aka  (Bka) 
annihilates  a  state  with  spin  a  and  momentum  k  located 
at  A  ( B )  site  [A  and  B  constitute  a  unit  cell  as  shown 
in  Fig.  2(c)].  Here,  ht  =  2J[cos{kxa)  +  cos {kya)]: 
bx  =  2Jso  sin(kxa)  and  by  =  —  2  Jso  sin(fcya);  tXiV  are 
Pauli  matrices  acting  on  A,  B  sublattices.  This  Hamil¬ 
tonian  can  be  transformed  into  a  block  diagonal  matrix, 
i.e.,  H  ->•  H'  =  (hz  +  ij)az  -  htcrzTz  +  rz{bxay  +  byax ), 
which  commutes  with  tz.  Note  that  we  have  neglected 
the  spin-independent  term  /io,  which  has  no  essential 
effects  on  the  physics. 

Similar  to  the  toy  model  in  Eq.  (1),  eigenvalues  of 
this  Hamiltonian  are  Eg±( k)  =  yjb±  —  y2  +  2ibz±'y  = 

yj A±{ k)eI0±/2,  where  A±(k)  =  ^ (b±  —  y2)2  +  46^±y2 

with  b±  =  b2x  +  b2  +  b2z±  and  bz±  =  ±ht  +  hz  (±  la¬ 
bel  two  particle  or  hole  bands  associated  with  the  sub¬ 
space  tz  =  - F  for  H')7  and  6±  are  defined  by  cos0±  = 
(6±  —  y2)/A±(k)  and  sin0±  =  26z±y/A±(k).  With¬ 
out  y,  energy  is  purely  real,  and  Weyl  points  emerge 
at  kxy0  =  (kxa,  kya,  kzaz)  =  [77  0,  —  2mTthz/ (h2k\z)]  or 
k w±  =  [0, 0, -2m7r(±4 J  +  hz)/(h2k2Lz)],  where  az  = 
Tt/kLz-  The  touching  point  is  fourfold  (doubly)  degen¬ 
erate  at  kwo  (kw±)-  When  y  >  0,  the  spectrum  be¬ 
comes  complex  and  it  is  purely  real  only  in  the  plane 
bz±  =  0.  A  touching  point  transforms  into  a  closed  line 
(i.e.,  Weyl  exceptional  ring)  at  which  particle  and  hole 
bands  coalesce  when  bz±  =  0  and  b2  +  b2  =  y2,  as  shown 
in  Fig.  3(a).  Around  k^0,  the  fourfold  degeneracy  of 
the  touching  point  is  broken,  and  there  emerge  two  Weyl 
exceptional  rings  that  are  not  degenerate  except  at  four 
points  with  |  sva.kxax\  =  y/(2y/2 Jso):  kx  =  ±fcy  —  77  and 
kz  =  — hz/a  [as  shown  in  Fig.  3(a)].  Around  k w±,  each 
Weyl  point  morphs  into  a  single  Weyl  exceptional  ring, 
which  can  be  approximated  by  k2  +  k2  =  y2/ (4  J|G)  and 
kz  =  [~hz  ±  J(4-  y2/(4J|0))]/a  when  y  <  2 JSo- 

Analogous  to  the  toy  model,  a  Weyl  exceptional  ring 
in  this  realistic  model  can  be  characterized  by  the  Chern 
number  defined  in  Eq.  (3),  i.e.,  evaluated  by  an  integral 
of  the  Berry  curvature  over  a  closed  surface  S  via  the 
Riemann  surface.  Around  kwoi  there  are  two  Weyl  ex¬ 
ceptional  rings  associated  with  two  branches  9± ,  and  the 
Chern  number  is  defined  for  each  band  with  Cg±  =  1 
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FIG.  3.  (Color  online)  (a)  Schematic  of  Weyl  exceptional  rings 
denoted  by  closed  red,  green  and  cyan  lines  for  the  system 
described  by  the  Hamiltonian  in  Eq.  (8).  The  dashed  box 
depicts  the  first  Brillouin  zone,  (b)  Real  (c)  imaginary  and 
(d)  absolute  values  of  the  eigenenergy  with  respect  to  kzaz 
for  ky  =  0  and  7  =  0.7 J  when  the  open  boundary  condition 
is  imposed  along  the  x  direction.  The  red  lines  are  the  surface 
states.  In  (d),  additional  surface  states  for  7  =  0,  0.35J,  0.86J 
are  plotted  as  blue,  green  and  yellow  lines  respectively.  Note 
that  only  the  parts  with  zero  absolute  energy  are  associated 
with  surface  states. 


[Cg±+2-K  =  —1)  when  S  encloses  one  ring.  Around 
kw±,  the  corresponding  band  contributes  Cg±  =  —  1 
(Ce±+2*  =  1).  Furthermore,  apart  from  the  Chern  num¬ 
ber,  we  can  characterize  the  ring  by  a  quantized  Berry 
phase  defined  in  Eq.  (6),  i.e.,  evaluated  along  a  closed 
trajectory  enclosing  the  Weyl  exceptional  ring  twice  for 
a  considered  band  with  a  ring. 

Another  intriguing  feature  of  Weyl  semimetals  is  the 
existence  of  a  Fermi  arc,  surface  states  that  connect  two 
Weyl  points  with  opposite  Chern  numbers  in  a  geometry 
with  edges.  When  7  >  0,  Weyl  points  develop  into  Weyl 
exceptional  rings,  and  one  may  wonder  on  the  existence 
of  surface  states  with  open  boundaries.  Here,  we  calcu¬ 
late  the  spectra  of  the  open  system  under  open  bound¬ 
aries  along  the  x  direction  and  plot  the  real,  imaginary 
and  absolute  parts  of  the  spectra  in  Fig.  3(b),  (c),  and 
(d),  respectively.  We  neglected  the  spin-independent  en¬ 
ergy  h2k2/2m  for  clarity.  Interestingly,  zero  energy  states 
emerge  for  both  real  and  imaginary  parts  of  the  spectra. 
Yet,  the  surface  states  (Fermi  arc)  are  only  associated 
with  those  of  zero  absolute  energy,  which  connect  the 
Weyl  exceptional  ring  at  the  center  (kp/0)  to  those  on  two 
sides  (kw±).  These  states  are  doubly  degenerate  eigen¬ 
vectors  (not  generalized  ones  [46,  47]),  one  (the  other) 
of  which  is  localized  on  the  left  (right)  surface.  Com¬ 
pared  to  the  surface  states  without  7,  their  range  along 
kz  decreases  with  respect  to  7  because  the  size  of  the 


rings  along  the  z  direction  grows  with  7.  Fig.  3(d)  shows 
the  shrinking  surface  states  for  7  =  0,  0.35J,  0.7 J  and 
0.86J.  As  7  becomes  sufficiently  large,  the  rings  around 
(kxax,kyay)  =  (0,0)  overlap  with  those  around  (77  0)  in 
the  kz  direction  and  surface  states  completely  disappear. 

To  measure  the  Weyl  exceptional  ring,  a  possible  ap¬ 
proach  is  to  probe  the  dynamics  of  atom  numbers  of  each 
spin  component  after  a  quench  [40] .  Initially,  if  we  only 
keep  the  spin-independent  optical  lattices  but  switch  off 
the  spin-dependent  ones  (contributing  to  the  spin-orbit 
coupling)  and  dissipation,  we  can  load  spin  up  atoms 
into  the  system  and  the  ground  state  is  ^(k  =  0,t  = 
0)  =  (  1  0  1  0  )  / x/2  since  the  Hamiltonian  reduces  to 
H  =  — htrx .  This  state  can  be  driven  to  a  state  with 
k  7^  0  by  accelerating  the  optical  lattices  or  by  apply¬ 
ing  an  external  gradient  force.  After  that,  the  spin-orbit 
coupling  and  dissipation  can  be  suddenly  turned  on.  So 
this  state  is  no  longer  the  eigenstate  of  the  system  and 
the  atom  numbers  will  change  with  time.  For  example,  if 
k  lies  in  the  bz_  =  0  plane,  the  normalized  atom  number 
for  spin  down  is  given  by 


ni  = 


b$ 


4|£e|2 


[£■ 

A— ± 


yX2lm(EQ)t/h 


—2  cos 


/  2R,  \ 


>  (9) 


where  nj.  =  7V^e7*/R  with  N±  being  the  atom  num¬ 
ber.  Outside  of  the  ring,  Im(Efi)  =  0  and  n j.  = 
(b2  +  b2)  sin2 (Egt/H)/ Eg  with  an  oscillation  period  of 
2nh/Eg  and  inside  the  ring  R e(Eg)  =  0  and  n j.  = 
{b2x  +  b2)[J2\=±  eA2Im('Be)t/?l  -  2]/(4|£e|2)  with  no  os¬ 
cillation.  The  existence  of  the  Weyl  exceptional  ring 
will  be  manifested  through  the  change  in  oscillation  pe¬ 
riods.  In  experiments,  one  may  choose  87Rb  (bosons) 
atoms  and  apply  blue-detuned  laser  beams  at  wavelength 
A  =  767  nm  [48]  to  generate  the  optical  lattices  with 
Weyl  points.  With  specific  experimental  settings,  our 
model  parameters  are  given  by  Jgo  =  0.5J  and  J  = 
0.058E(r,  where  the  recoil  energy  is  Er/H  =  hkji/2m  = 
2tt  x  3.9  kHz  with  kn  =  2tt/X  and  A  being  the  wavelength 
of  laser  beams.  The  decay  strength  7  can  be  experimen¬ 
tally  tuned  by  controlling  the  intensity  of  the  resonant 
optical  beam. 

In  summary,  we  have  discovered  a  Weyl  exceptional 
ring  in  a  dissipative  system  of  Weyl  points  with  particle 
gain  and  loss.  Such  a  ring  is  characterized  by  both  a 
quantized  Chern  number  and  a  quantized  Berry  phase, 
which  are  defined  via  the  Riemann  surface.  We  further 
propose  an  experimental  scheme  in  cold  atoms  to  realize 
the  Weyl  exceptional  ring,  which  paves  the  way  for  future 
experimental  investigation  of  such  a  ring  and  its  unusual 
topological  properties. 
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SUPPLEMENTAL  MATERIAL 


In  the  supplemental  material,  we  will  derive  the  semiclassical  equation  for  a  wave  packet  in  an  optical  lattice  with 
a  non-Hermitian  term.  The  wave  packet  can  be  constructed  from  Bloch  wave  functions: 

\&)  =  J  dka{k,t)\<f>n(k)),  (SI) 

where  |<^„(k))  is  a  right  Bloch  wave  function  of  He,  the  Hamiltonian  describing  the  system,  i.e.,  He |</>„(k))  = 
En{k)\(j)n{k))  with  |<^„(k))  =  elk  r|u„(k)),  (u„(k)|u„(k))  =  1  and  n  being  the  band  index;  a(k,i)  =  |a(k, t)|e*7^k'^ 
with  the  amplitude  |a(k,  t)\  taking  the  Gaussian  form  centered  at  kc  =  J  dk|a'(k,  f)|2k  with  |a,' (k,  f)  |2  = 
|a(k,  t)\2 /  f  dk'|a(k',  t)|2  and  the  phase  j(k,t)  incorporating  a  dynamical  and  Berry  phase.  Since  we  only  consider 
the  case  that  a  single  band  is  occupied  by  the  wave  packet,  we  will  neglect  the  subscript  n  hereafter. 
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Using  Eq.  (SI),  the  location  of  the  wave  packet  is  given  by 

(52) 

(53) 


rc=  (<f>|r|$)/($|$) 

=  —J  dk\ci'(k,  t)|29k7(k,  t)  +  j  dfc|a'(k,  t)|2(u(k)|z<9k|ii(k)) 

~  -<9kc7 (kc,t)  +  A(kc), 


where  A(kc)  =  z(zz(kc)|<9kc  |zt(kc))  and  7(kc,t)  =  -  ftg  dt'Re  Ek^tt)  -  A(kc(t'))  •  y|f  with  A(kc)  = 

i(u(kc)|9k<,|u(kc)),  where  the  first  and  second  terms  correspond  to  the  dynamical  and  Berry  phases,  respectively. 
Using  Eq.  (S3),  we  can  obtain  the  semiclassical  equation  in  terms  of  a  velocity 


z(to)  =  lim5t^0- 


„(to  +  St)  -  rc(t0) 
St 


=  dkcE( kc)  -  kc  x  fl(kc) 


(S4) 


with  E( kc)  =  Re[U(kc)]  +  A(kc)  •  kc. 


